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Despite the spontaneity of some in vitro protein folding reactions, native folding in vivo often 
requires the participation of barrel-shaped multimeric complexes known as chaperonins. Although 
it has long been known that chaperonin substrates fold upon sequestration inside the chaperonin 
barrel, the precise mechanism by which confinement within this space facilitates folding remains 
unknown. In order to address this question, we focused on the solvent of the folding reaction and 
investigated the effect of confinement within a chaperonin-like cavity on the configurational free 
energy of water. Using a simple theory for the thermodynamics of density and hydrogen bond 
order fluctuations in liquid water, we calculated solvent free energies for cavities corresponding to 
the different conformational states in the ATP-driven catalytic cycle of the prokaryotic chaperonin 
GroEL. Our findings suggest that one function of chaperonins may be to trap unfolded proteins 
and subsequently expose them to a micro-environment in which the hydrophobic effect, a crucial 
thermodynamic driving force for folding, is enhanced. 



INTRODUCTION 

Efficient production of natively folded protein is an 
absolute requirement for cell survival. Although the 
amino acid sequence of a small polypeptide generally suf- 
fices to specify the chain's native structure under condi- 
tions favorable to spontaneous folding in vitro, the pro- 
tein folding problem in vivo is more challenging in several 
respects. To begin with, the cell must somehow protect 
the delicate homeostasis of its proteome against various 
environmental stresses that wreak havoc with the ther- 
modynamics and kinetics of folding. Moreover, folding 
reactions that take place in the highly crowded intracellu- 
lar milieu run the risk of never proceeding to completion 
because of the tendency of partially unfolded polypep- 
tides to aggregate non-specifically and irreversibly at 
high concentrations. And even under more dilute con- 
ditions, a large protein that circuitously explores many 
non-native kinetic traps along the way to its native state 
may fold too slowly to be of any use in a biological con- 
text. 

Proteins known as molecular chaperones enable the 
cell to overcome these obstacles [TJ [5]. Chaperones are 
characterized by their ability to bind the exposed hy- 
drophobic moieties of non-native polypeptides and guide 
these wayward chains back onto productive folding path- 
ways. Perhaps the most intriguing type of chaperone is 
the chaperonin, a barrel-shaped multimeric complex that 
engulfs and releases its substrates in an ATP-dependent 
manner [2j |3j H]- Chaperonins sub-divide into groups 
I (eubacterial) and II (archaebacterial/eukaryotic), and 
the group I tetradecamer Hsp60/GroEL from E. coli is 
the most extensively studied model of chaperonin action. 
The open conformation of GroEL uses its largely hy- 



drophobic inner surface to snare an unfolded substrate. 
Subsequent binding of ATP and the GroES co-factor in- 
duces a conformational change in GroEL, after which the 
substrate may attempt to fold inside a closed complex 
whose interior cavity surface displays many hydrophilic 
residues. Once ATP hydrolysis takes place, GroES disso- 
ciates from the complex and the barrel re-opens, releasing 
the substrate into the surroundings and completing the 
reaction cycle. 

For some of its substrates, GroEL mainly provides 
an infinite-dilution environment in which a protein may 
fold without risk of aggregation. For others, however, it 
has been shown that confinement within the chaperonin 
barrel is accompanied by a marked increase in the rate 
of folding even under conditions that preclude aggrega- 
tion [5j. The origins of this chaperonin foldase activity 
are still largely unclear, but recent experiments have be- 
gun to shed some light on the question [6J. Tang et. 
al. assayed the activity of a wide range of GroEL mu- 
tants and found that the acceleration of folding achieved 
by a chaperonin mutant correlated with the net charge 
and hydrophilicity of its interior cavity surface. This in- 
triguing result suggests that the hydrophilic walls of the 
cavity somehow reshape the folding energy landscape of a 
confined protein to favor progression to the native state 
[1]. At present, the mechanism for this reshaping re- 
mains mysterious, and an explanation for how it might 
come about would undoubtedly do much to improve our 
understanding of chaperonin action. 

This study aims to provide such an explanation by 
considering the role of water in folding under confine- 
ment. Folding reactions take place in aqueous solvent, 
and it is the strong tendency of water molecules to hydro- 
gen bond with each other that generates the "hydropho- 
bic effect" |7j, a thermodynamic force that helps drive 



proteins to fold by causing non-polar amino acid residues 
to reduce their solvent-exposed surface area [S, • Hydro- 
gen bond quality is sensitive to the relative orientations of 
molecules, and numerous studies have documented sub- 
stantial changes in the behavior of water once the liquid 
is confined in a small enough space that the organization 
of the hydrogen bond network becomes significantly con- 
strained [9l [ini [11] . It is therefore reasonable to hypoth- 
esize that the high degree of confinement experienced by 
water that participates in folding reactions inside GroEL 
may substantially alter the solvent free energy, and there- 
fore the folding landscape, relative to what it would be 
in the cytosol. 

Here, we investigate the capacity of a nanopore to 
modulate the free energy of water confined within it by 
undergoing shifts in surface chemistry that mimic the 
conformational changes in GroEL 's catalytic cycle. Using 
a phenomenological theory for the free energy of fluctua- 
tions in density and hydrogen bond order in bulk water, 
we compute the effect of confinement in hydrophilic and 
hydrophobic cavities on the "folding landscape" of a sim- 
ple model protein. Our findings suggest that GroEL may 
be designed to preferentially bind unfolded substrates 
and allow them to fold in a micro-environment in which 
the hydrophobic effect is enhanced above cytosolic levels. 

MODEL 

In order to examine the possible role of GroEL in re- 
modeling the solvent of the folding reaction, we sought to 
develop an analytical framework for describing the ther- 
modynamics of confined water. Microscopic theories of 
molecular liquids are notoriously difficult, especially for 
a molecule with water's rare tendency to form extensive, 
tetrahedrally-coordinated networks of hydrogen bonds. 
For this reason, we chose to avoid microscopies in favor 
a phenomenological approach. Following past work |12j 
that attempted to explain the anomalous bulk properties 
of water using a Landau-like statistical field theory, we 
posited that the effective free energy of a fluctuation in 
the configuration of bulk water at a given temperature T 
and chemical potential ji depends entirely on two order- 
parameters: the local density of the fluid p(r) and the 
local quality of the hydrogen bond network s(r). Defin- 
ing 5p and 5s to be the respective deviations of these 
quantities from their average values in the bulk, we can 
write the free energy of water for small configurational 
fluctuations in a volume V as 

m,5s^ + TUpSp^ - e 5s 5p] + 0{5^) (1) 

The coefficients multiplying each term in the integral 
should, in principle, be functions of /x and T, although we 



would expect the coupling e between the density and hy- 
drogen bond order to roughly equal the energy of a hydro- 
gen bond. The derivative terms simply ensure that both 
order parameters not vary too rapidly over space, since 
they are only well-defined quantities on length scales 
larger than a water molecule. The squared terms ensure 
that larger deviations from bulk values of the order pa- 
rameter carry higher free energetic costs, while the cross- 
term proportional to 5p 5s dictates that the formation of 
more hydrogen bonds lowers the energy of the system. 

It is important to note here that because we did not 
include cubic and quartic terms in our free energy func- 
tional, we would be unable to use it to capture the trade- 
off between fluid density and hydrogen bond quality that 
leads to the anomalous decrease in the density of water 
as it cools to its freezing point. While such effects may, 
in fact, be important to adequately describing conflne- 
ment effects in water, our initial analysis assumes for the 
sake of simplicity that the deviations from bulk behavior 
involved in the phenomenon in question are small enough 
that the quadratic functional written above will suffice. 

Thus far, we have presented a simple theory of or- 
der fluctuations in bulk water. Before we can apply this 
theory to the study of confined water, however, we must 
augment it by introducing an additional contribution to 
the free energy functional that comes from the interac- 
tion between the liquid and the surface that encloses it. 
Thus, let our system consist of water confined within a 
shell S of thickness w, where w is comparable to the 
size of a single water molecule. If S is perfectly non- 
polar, then it cannot participate in hydrogen bonds and 
we must require that the contribution to the free energy 
due to hydrogen bonding that would be given by e p s in 
the bulk must be zero in the shell. We may further in- 
troduce some potential e u(5p, 5s) that accounts for any 
additional effects from the surface chemistry of S. Thus, 
we obtain a new term in our functional: 

J-s — we dS [X p s + u{p, s)] (2) 
Js 

Here, < A < 1 is a parameter that reflects the degree of 
hydrophobicity of the surface, with A = 1 corresponding 
to a completely non-polar surface. 

To find the equilibrium density and hydrogen bond 
order, we minimize the total free energy J^t = + 
with respect to the order parameter fields by using the 
calculus of variations. The equilibrium fields 6peq and 
5seq will obey a system of Helmholtz equations that may 
be de-coupled by a change of basis. In the cases of spheri- 
cally and planar symmetric systems, it is straightforward 
from there to obtain closed-form solutions whose inte- 
gration constants must be set by plugging back into the 
functional and minimizing J-'x- For arbitrary cavity ge- 
ometries and surface chemistries, however, it is necessary 
to seek a numerical solution. We therefore restrict our 



attention here to the analytically tractable cases, which 
are somewhat easier to deal with and nevertheless suffice 
as tools for examing the issues of interest in this study. 



RESULTS 

We modeled GroEL as a spherical pore of radius Rq 
and the chaperonin substrate protein as a smaller sphere 
of radius Rp. The chaperonin barrel was assumed to have 
two different states that reflected the open and closed 
conformations identified in structural studies of GroEL 
[IB] . In the open state, the inner surface of the chaper- 
onin was treated as a largely non-polar surface, meaning 
Aopen = 1 and Uopen{Sp,5s) = 0. In the closed state, 
the chaperonin cavity was treated as being highly hy- 
drophilic because of its richness in charged and polar 
residues that participate in hydrogen bonds with water 
molecules. Thus, we set Xdosed — and Udosed — —u 5p, 
where u > accounts for the net attraction of favorably- 
oriented molecular electric dipole moments to the net 
charge of the cavity surface. The protein was also as- 
sumed to have two states: native and unfolded. The 
native conformation was assumed to have a hydrophilic 
surface chemistry described by Xfoided = Ufoided = 0, 
and the unfolded state was assumed to be partially hy- 
drophobic, with Uunfolded = and Xunfolded > 0. 

The solutions to the resulting Helmholtz equa- 
tions were linear combinations of terms of the form 
aexp[±cr]/r, with the coefficients a determined through 
minimization of the total free energy J-'t, and the cor- 
relation length scales 1/c set by the parameters of the 
theory. In order to plot these solutions, it was neces- 
sary to choose specific values for the parameters in the 
free energy functional. Since the questions we sought to 
answer were qualitative in nature, the precise parame- 
ter set used for the study was somewhat arbitrary. In- 
deed, we were unable to observe any qualitative change 
in the results from varying the coefficients so long as it 
was assured that bulk density and hydrogen bond order 
remained a free energy minimum. The constraints we 
ultimately made sure to satisfy were that the shell width 
w and the characterstic correlation length scales for the 
order parameters 1/c were comparable to the size of a 
few water molecules and that the protein and chaperonin 
cavity were sized correspondingly. 

Figure [l] displays the equilibrium liquid density and 
hydrogen bond order profiles for the different states of 
GroEL and its substrate. For the case of a non-native 
substrate inside an open GroEL complex, both surfaces 
confining the solvent are hydrophobic (Fig. [l] top). As 
a result of the loss of hydrogen bonding at these sur- 
faces, water withdraws from them and from the cavity 
as a whole, driving the average density below bulk lev- 
els. Moreover, it is apparent from a calculation of the 
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FIG. 1: Equilibrium order parameters Speq and Sseq inside 
spherical shells were calculated using nip = 20, nis = 10, 
Dp = 20, Ds = 30, e — 5, po — So = 1, w = .5, u — 1, Rp — 3, 
and Ra = 7. All data reported in this study use the same 
parameters for the free energy functional. Left column: Liq- 
uid density (darker shading corresponding to lesser density) 
is plotted in the space between the surface of the protein (red 
unfolded, green folded) and the cavity (red hydrophobic, blue 
hydrophilic) wall for an unfolded protein inside an "open" hy- 
drophobic cavity (top), an unfolded protein inside a "closed" 
hydrophilic cavity (middle), and a folded protein in a closed 
cavity (bottom). Gold arrows indicate the direction of the 
solvent-mediated force between the two surfaces (See Fig. [2]|. 
Right column: Liquid density (green curves) and hydrogen 
bond order (orange curves) are plotted as a fraction of their 
bulk values. At a hydrophobic surface there is a loss of hy- 
drogen bonding and a depletion of liquid density. In contrast, 
at a highly hydrophilic surface there is an elevation in liquid 
density and a greater amount of hydrogen bonding. 

free energy of solvent trapped between two planar sur- 
faces of the same type that the force between the non- 
native substrate and the chaperonin wall is attractive 
(Fig. [2] red curve), and would favor binding. Once 
the chaperonin forms a closed complex and 89890under- 
goes its conformational change, however, the projection 
of charged and polar side-chains into the barrel's interior 
causes the solvent order and resulting thermodynamic 
forces to change. With the inner surface of the barrel 
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FIG. 2: Solvent free energy in arbitrary units is plotted as 
a function of separation L for pairs of horizontal plates. Dis- 
tance is measured in units of the shorter correlation length 
of the Landau theory, roughly equivalent to the size of a few 
water molecules. Between two hydrophobic plates (red curve) 
the solvent mediates an attractive force that grows stronger 
with proximity. Between a hydrophobic plate and a highly hy- 
drophilic plate (blue curve), there is a weaker repulsive force. 
Between two hydrophilic plates the force is essentially non- 
existent (green curve). 



now highly hydrophilic, the solvent density is elevated 
near the chaperonin wall and depressed surrounding the 
unfolded protein (Fig. [l] middle) . The result is a repul- 
sive force between the surfaces that pushes the substrate 
towards the center of the cavity (Fig. |2] blue curve). 
Upon folding, the substrate surface becomes more hy- 
drophilic, allowing the density of the nearby solvent to 
relax to bulk levels (Fig. [l] bottom) with the result that 
the force between chaperonin and substrate becomes very 
weak (Fig. [2] green curve). 

With equilibrium density and bond order profiles in 
hand, we were able to compute the solvent free energy 
changes for folding inside the different conformers of 
GroEL (Fig. |3]) . Within the framework of the model con- 
sidered here, folding without the involvement of a chaper- 
onin simply consists of a gradual reduction in hydropho- 
bicity of the surface of a sphere of radius Rp surrounded 
by aqueous solvent. In contrast, chaperonin- assisted fold- 
ing divides into several stages. First, the hydrophobic 
interior of the barrel in its open conformation decreases 
its free energy by binding to hydrophobic moeities ex- 
posed by a partially unfolded substrate. Once the non- 
native protein is inside the cavity, the solvent free energy 
change associated with folding is made less negative by 
the surrounding hydrophobic walls, and the chaperonin 
therefore provides a thermodynamic drive towards un- 
folding that increases as the radius of confinement be- 
comes smaller (Fig. [Sj red curve). This effect comes 
about because the hydrophobic cavity walls deplete the 
amount of water in the cavity, making the cost of being 
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FIG. 3: The change in the folding free energy A AG/ is de- 
fined as the difference between the free energy of folding un- 
der confinement AG} and the free energy of folding in bulk 

solution AG^"' . The fractional change in this solvent free en- 
ergy of folding is plotted as a function of confinement radius 
Ai? = Rg — Rp for _Rp = 3 less the width w of the inter- 
action shell at each surface. Length is measured in units of 
the shorter correlation length of the theory of solvent fluctua- 
tions described in equationjl] When the protein is confined in 
the largely hydrophobic open GroEL cavity (red curve), the 
folded state is destabilized moreso as the degree of confine- 
ment increases. However, when confined in the closed, highly 
hydrophilic chaperonin cavity (blue curve), stabilization of 
the folded state increases with the degree of confinement. 

unable to form hydrogen bonds at the protein surface 
lower than in bulk solution. After GroEL binds to GroES 
and ATP and undergoes its conformational change, how- 
ever, the reverse is true: the solvent portion of the free 
energy of folding inside the cavity becomes more nega- 
tive than that of folding in bulk solution, with the mag- 
nitude of the effect again increasing as the confinement 
becomes more severe (Fig. [s] blue curve). This dra- 
matic shift takes place because the charged inner surface 
of the chaperonin has elevated the density of the water 
in the cavity, thereby raising the number of molecules 
who have difficulty finding a hydrogen bonding partner 
when hydrophobic substrate surface area is presented to 
the solvent. Thus, our model predicts that the confor- 
mational change in GroEL brings about a remodeling of 
the confined solvent that enhances the hydrophobic effect 
above its strength in bulk solvent and helps to drive the 
folding reaction to completion. 



DISCUSSION 

Past attempts to explain the foldase activity of chap- 
eronins have pointed to a diverse list of possible causes. 
While some researchers have focused on evidence that 
GroEL may assist in the unfolding of kinetically trapped 




FIG. 4: A model of chaperonin action. Hydrophobic forces 
cause an unfolded (red) protein to bind to the wall of the 
open (red) GroEL barrel {a ~> b). Upon formation of a closed 
(blue) complex containing GroES and ATP, rearrangements 
in the barrel present a hydrophilic surface to the interior that 
repels the substrate into the center of the cavity {b ^ c). Be- 
cause of this repulsion, the free energy AG*f of folding to the 
native (green) conformation inside the cavity (c d) is more 
negative than the folding free energy AGj?' out in bulk solu- 
tion (a e). The hydrolysis of ATP therefore drives a local 
enhancement of the hydrophobic effect inside the chaperonin. 

substrates [Ml [15], others see a role for the barrel-shaped 
cavity in accelerating folding by reducing the conforma- 
tional entropy of a protein's unfolded state [TB], while still 
others propose that interactions between the protein and 
side-chains on the cavity wall might help to smooth out 
some of the pitfalls in the protein's free energy landscape 
[171 118] . In this work, we examined a novel dimension 
of the chaperonin puzzle by studying the thermodynam- 
ics of the aqueous solvent in which the folding reaction 
takes place. By developing a simplified theory of order 
fluctuations in confined water, we were able compute the 
free energy of water in cavities that mimicked the salient 
features of a protein folding inside GroEL. Our findings 
underlined the potential importance of solvent effects in 
a complete description of how GroEL promotes folding. 

Figure 4 summarizes the model of chaperonin ac- 
tion that emerges from our study. A non-native protein 
that exposes hydrophobic groups to the bulk solvent will 
be engulfed by an open GroEL complex because of hy- 
drophobic forces mediated by the solvent. At this point, 
the substrate becomes bound to the wall of the chaper- 
onin. Once GroES and ATP bind, however, the solvent 
travels uphill in free energy as a conformational change in 
GroEL presents charged residues to the barrel's interior. 
As a result, the substrate is forced away from the walls 
into the center of the cavity, a prediction that is consis- 
tent with experimentally measured increases in substrate 
fluorescence anisotropy following closing of the complex 
[6]. The subsequent folding reaction takes place inside 
a charged, hydrophilic chamber, and is therefore accom- 
panied by a more drastic decrease in the free energy of 
the surrounding solvent. In other words, the interior of 



GroEL may provide a micro-environment in which the 
hydrophobic effect is substantially stronger than in bulk 
solution. 

This finding leads us to propose a novel mechanism 
for the acceleration of folding by chaperonins. The hy- 
drophobic effect is known to be a crucial contributor 
to the thermodynamic stability of natively folded pro- 
teins. However, several studies have also documented 
the importance of hydrophobic forces to folding kinet- 
ics [m [2Ql [21], both in terms of bringing about nucle- 
ation through hydrophobic collapse [5D] and through sta- 
bilization of the folding transition state [T9j. It is there- 
fore plausible that the stronger hydrophobic effects ex- 
perienced by substrates inside closed GroEL complexes 
would hasten progression towards the native state. Thus, 
through its interaction with the confined solvent of the 
folding reaction, the chaperonin could reshape the free 
energy landscape in order to promote more rapid fold- 
ing. 

This hydrophobic enhancement model of chaperonin 
action provides us with a ready means to explain the 
puzzling relationship between foldase activity and cav- 
ity hydrophilicity reported in recent experiments. Our 
model assumes that charged groups on the interior of 
GroEL draw water into the barrel chamber, and that 
this elevation of solvent density leads directly to an en- 
hancement of hydrophobic effects. The model therefore 
predicts that replacement of charged amino acid residues 
with neutral ones inside the GroEL complex should sig- 
nificantly reduce the chaperonin 's foldase activity, just 
as was observed in the recent mutational studies of Tang 
and co-workers [Bj. Interestingly, the same study also 
reported a negative correlation between foldase activity 
and the volume of the GroEL folding chamber. Although 
this result was interpreted to be due to the reduction in 
conformational entropy of the unfolded protein that con- 
finement causes, it is worth noting that since hydrophobic 
enhancement rises with the degree of confinement (Fig. 
3), our model provides an alternative explanation for the 
outcome of those experiments as well. 

Additional evidence pointing to charged surface 
residues as an important determinant of GroEL activ- 
ity comes from the work of Wang et. al. [22 , in which 
GroEL-GroES was optimized through directed evolution 
to fold green fluorescent protein (GFP) more rapidly. 
GFP is a small protein with an exceptionally well-buried 
hydrophobic core that would be likely to form more 
rapidly in response to enhancement of hydrophobic ef- 
fects. The authors found that a key mutation for accel- 
erating the folding of GFP was the replacement of an 
aromatic tyrosine on the interior surface of GroES with 
either arginine or histidine, both charged residues. These 
results are consistent with our proposal that charged cav- 
ity residues accelerate folding by strengthening the hy- 
drophobic effect. 



The existing evidence consistent with the mechanism 
of chaperonin foldase activity proposed here is substan- 
tial, but too indirect to be compelling by itself. For- 
tunately, our model's most endearing quality is that it 
makes clear predictions that can easily be tested in sil- 
ica and in the laboratory. At the most basic level, we 
would expect an all atom simulation of water and hy- 
drophobic solutes inside a sufficiently hydrophilic cavity 
to show that hydrophobic forces between the solutes are 
enhanced by the confinement. We would also anticipate 
that a full-scale simulation of the closed GroE complex 
would show that water density inside the cavity was el- 
evated (a prediction already confirmed by unpublished 
results), and that the degree of the corresponding eleva- 
tion in simulations of previously characterized charged- 
to-neutral mutants would correlate roughly with their 
measured foldase activity. Finally, and most conclusively, 
we would predict that experiments that used Forster res- 
onance energy transfer (FRET) or other techniques to 
measure the strength of hydrophobic forces between so- 
lutes inside GroEL would demonstrate that closure of the 
barrel through the binding of GroES and ATP is nec- 
essary and suSicient to enhance the hydrophobic efi^ect 
inside the cavity. 
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